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Abstract 

We recall the physical features of the parton distributions in the quantum 
statistical approach of the nucleon, which allows to describe simultaneously, 
unpolarized and polarized Deep Inelastic Scattering data. Some predictions 
from a next-to-leading order QCD analysis are compared to recent exper- 
imental results and we stress the importance of some tests in the high-x 
region, to confirm the validity of this approach. 
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Let us first recall some of the basic ingredients for building up the parton 
distribution functions (PDF) in the statistical approach, as oppose to the 
standard polynomial type parametrizations, based on Regge theory at low x 
and counting rules at large x. The fermion distributions are given by the sum 
of two terms [Ij, the first one, a quasi Fermi-Dirac function and the second 
one, a flavor and helicity independent diffractive contribution equal for light 
quarks. So we have, at the input energy scale Ql = 4GeV^, 

xq {X, Q,) - ^^^^^^ _ ^ ^ + ^^p^^/-^ ^ ^ , (Ij 
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xq (X, go) - ^^p^^^ ^ ^-Hy-^ ^ 1 + exp(x/x) + 1 ' 

Notice the change of sign of the potentials and helicity for the antiquarks. 
The parameter x plays the role of a universal temperature and X^^ are the 
two thermodynamical potentials of the quark g, with helicity h = ±. It 
is important to remark that the diffractive contribution occurs only in the 
unpolarized distributions q{x) = q+{x) +g_(x) and it is absent in the valence 
qv{x) = q{x) — q{x) and in the helicity distributions Ag_(x) = q+{x) — 
(similarly for antiquarks). The eight free parameterqj in Eqs. ( TTPj) were 
determined at the input scale from the comparison with a selected set of 
very precise unpolarized and polarized Deep Inelastic Scattering (DIS) data 
[1]. They have the following values 

X = 0.09907, b = 0.40962, b = -0.25347, A = 0.08318, (3) 

X(+ = 0.46128, = 0.29766, X^" = 0.30174, X^^ = 0.22775 . (4) 
For the gluons we consider the black-body inspired expression 

xG{x,Ql) = -, (5) 

exp(x/x) — 1 

a quasi Bose-Einstein function, with bo = 0.90, the only free parameter 
since Ac = 20.53 is determined by the momentum sum rule. We also assume 
that, at the input energy scale, the polarized gluon distribution vanishes, 
so xAG{x, Ql) = 0. For the strange quark distributions, the simple choice 
made in Ref. pTj was greatly improved in Ref. [2]. More recently, new tests 
against experimental (unpolarized and polarized) data turned out to be very 
satisfactory, in particular in hadronic collisions, as reported in Refs. [H H] . 

For illustration, we will just give one recent result, directly related to 
the determination of the quark distributions from unpolarized DIS. We dis- 
play on Fig. l{Left), the resulting unpolarized statistical PDF versus x 
at Q2=io GeV^ where xUy is the M-quark valence, xo?^ the c?-quark va- 
lence, with their characteristic maximum around x = 0.3, xG the gluon 



A = 1.74938 and A — 1.90801 are fixed by the following normalization conditions 
u — u = 2, d — d — I. 

•^In Ref. [Tj we were assuming that, for very small x, xG{x, Qq) has the same behavior 
as xq{x, Qq), so we took bo = ^ + b. However this choice leads to a too much rapid rise 
of the gluon distribution, compared to its recent determination from HERA data, which 
requires be — 0.90. 
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Figure 1: Left : BBS predictions for various statistical unpolarized parton 
distributions versus x at Q"^ = lOGeV^. Right : Parton distributions at 
= lOGeV^, as determined by the HIPDF fit, with different uncertainties 
(Taken from Ref. [5]). 



and xS stands for twice the total antiquark contributions, i.e. xS{x) = 
2x{U{x) + d{x) + s{x)) + c(x)). Note that xG and xS are downscaled by a 
factor 0.05. They can be compared with the parton distributions as deter- 
mined by the HIPDF 2009 QCD NLO fit, shown also in Fig. l{Right), and 
the agreement is rather good. The results are based on recent ep collider data 
from HERA, combined with previously published data and the accuracy is 
typically in the range of 1.3 - 2%. 

Another interesting point concerns the behavior of the ratio d{x)/u{x), 
which depends on the mathematical properties of the ratio of two Fermi-Dirac 
factors, outside the region dominated by the diffractive contribution. So for 
a; > 0.1, this ratio is expected to decrease faster for Xq^—x < x < Xq^+x and 
then above, for x > 0.6 it flattens out. This change of slope is clearly visible 
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in Fig. |2] (Left), with a very little dependence. Note that our prediction 
for the large x behavior, differs from most of the current literature, namely 
d(x)/u{x) — )■ for X — >• 1, but we find d{x)/u{x) — )■ 0.16 near the value 
1/5, a prediction originally formulated in Ref. [6]. This is a very challenging 
question, since the very high-a; region remains poorly known. 




Figure 2: Left : The ratio d{x)/u{x) as function of x for = 4GeV^ (solid 
line) and = lOOGeV^ (dashed-dotted line). Right : Comparison of the 
data on d/u{x,Q'^) from E866/NuSea at = 54GeV^ [7], with the predic- 
tion of the statistical model (solid curve) and the set 1 of the parametrization 
proposed in Ref. [8] (dashed curve). 

To continue our tests of the unpolarized parton distributions, we must 
come back to the important question of the flavor asymmetry of the light 
antiquarks. Our determination of u{x, Q"^) and d{x, Q^) is perfectly consistent 
with the violation of the Gottfried sum rule, for which we found the value 
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Figure 3: The interference term xF^^ extracted in e^p collisions at HERA. 
Data from [TT] compared to the BBS prediction.. 

Ig = 0.2493 for = 4GeV^ Nevertheless there remains an open problem 
with the X distribution of the ratio d/u for x > 0.2. According to the Pauli 
principle, this ratio is expected to remain above 1 for any value of x. However, 
the E866/NuSea Collaboration [7] has released the final results corresponding 
to the analysis of their full data set of Drell-Yan yields from an 800 GeV/c 
proton beam on hydrogen and deuterium targets and they obtain the ratio, 
for = 54GeV^, d/u shown in Fig. |2] (Right). Although the errors are 
rather large in the high-x region, the statistical approach disagrees with the 
trend of the data. Clearly by increasing the number of free parameters, it 
is possible to build up a scenario which leads to the drop off of this ratio 
for X > 0.2. For example this was achieved in Ref. [S], as shown by the 
dashed curve in Fig. [2] (Right). There is no such freedom in the statistical 
approach, since quark and antiquark distributions are strongly related. On 
the experimental side, there are now new opportunities for extending the d/u 
measurement to larger x up to x = 0.7, with the upcoming E906 experiment 
at the 120 GeV Main Injector at Fermilab |H] and a proposed experiment at 
the new 30-50 GeV proton accelerator at J- PARC [10]. 
One can also test the behavior of the interference term between the photon 
and the Z exchanges, which can be isolated in neutral current e^p collisions 
at high Q"^. We have to a good approximation, if sea quarks are ignored. 
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xFf = f (2m„ + d^) and the comparison between data and prediction is 
displayed in Fig. [3l Here again, we note the remarkable agreement in the 
high- a; region. 
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Figure 4: Left : Quark and antiquark helicity distributions as a function of 
X for = 3GeV^. Data from COMPASS [I2]. The curves are predictions 
from the statistical approach. Right : Ratios (Am + Am) /{u + u) and [Ad + 
Ad)/{d + d) function of x. 
and a JLab Hall A experiment 
statistical approach. 



Data from Hermes for = 2.5GeV^ [13] 
The curves are predictions from the 



Analogous considerations can be made for the corresponding helicity dis- 
tributions, whose most recent determinations are shown in Fig. H] (Le/f). By 
using a similar argument as above, the ratio Au{x)/u{x) is predicted to have 
a rather fast increase in the x range {X^^ — x, Xq^ + x) and a smoother be- 
haviour above, while Ad{x) / d{x), which is negative, has a fast decrease in the 
X range {X^^ — x, X^^ + x) and a smooth one above. This is exactly the trends 
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displayed in Fig.\^{Right) and our predictions are in perfect agreement with 
the accurate high-a; data. We note the behavior near x = 1, another typical 
property of the statistical approach, is also at variance with predictions of the 
current literature. The fact that Au{x) is more concentrated in the higher x 
region than Ad{x), accounts for the change of sign of gi{x), which becomes 
positive for x > 0.5, as first observed at Jefferson Lab [H] . 

Concerning the light antiquark helicity distributions, the statistical ap- 
proach imposes a strong relationship to the corresponding quark helicity 
distributions. In particular, it predicts Au{x) > and Ad{x) < 0, with 
almost the same magnitude, in contrast with the simplifying assumption 
Au{x) = Ad{x), often adopted in the literature. The COMPASS experiment 
at CERN has measured the valence quark helicity distributions, defined as 
Aqy{x) = Aq{x) — Aq{x). These recent results displayed in Fig. [5] are com- 
pared to our prediction and the agreement is best in the high-x region. 
The data give Au{x) + Ad{x) ~ 0, which implies either small or opposite 
values for Au{x) and Ad{x). Indeed Au{x) > and Ad{x) < 0, predicted by 
the statistical approach [T] (see Fig. H] (Left), lead to a non negligible positive 
contribution of the sea to the Bjorken sum rule, an interesting consequence. 

We now turn to another important aspect of the statistical PDF and very 
briefly discuss a new version of the extension to the transverse momentum 
dependence (TMD). In Eqs. (HP the multiplicative factors Xg^^ and (X^g^)-^ 
in the numerators of the non-diffractive parts of g's and g's distributions, im- 
ply a modification of the quantum statistical form, we were led to propose in 
order to agree with experimental data. The presence of these multiplicative 
factors was justified in our earlier attempt to generate the TMD [16], but it 
was not properly done and a considerable improvement was achieved recently 
[17J. We have introduced some thermodynamical potentials Y^^, associated 
to the quark transverse momentum fcy, and related to Xq^ by the simple 
relation ln(l + exp[Fo^]) = kXl}g. We were led to choose k=3.05 and this 
method involves another parameter /i^, which plays the role of the tempera- 
ture for the transverse degrees of freedom and whose value was determined 
by the transverse energy sum rule. We have calculated the pt dependence of 
semiinclusive DIS double longitudinal-spin asymmetries, taking into account 
the effects of the Melosh-Wigner rotation, for vr^ production, which were 
compared to recent experimental data from CLAS at JLab. 
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Figure 5: The valence quark helicity distributions versus x and evolved 
at Q"^ = lOGeV^. The solid curve is the BBS prediction of the statistical 
approach and the data come from Ref. jl5] . 

A new set of PDF was constructed in the framework of a statistical ap- 
proach of the nucleon. All unpolarized and polarized distributions depend 
upon nine free parameters for light quarks and gluon, with some physical 
meaning. New tests against experimental (unpolarized and polarized) data 
on DIS, Semi- inclusive DIS and also hadronic processes, are very satisfac- 
tory. It has a good predictive power, but some special features remain to be 
verified, specially in the high-x region, a serious challenge for the future. 
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